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ABSTRACT 

Aims. In this paper we use a well-controlled spectroscopic sample of galaxies at l<z<3 drawn from the Galaxy Mass Assembly ultra¬ 
deep Spectroscopic Survey (GMASS) to study different star formation rate (SFR) estimators. In particular, we use infrared (IR) data to 
derive empirical calibrations to correct ultraviolet (UV) and [Oil]A3727 luminosities for dust extinction and dust-corrected estimates 
of SFR. 

Methods. We selected 286 star-forming galaxies with spectroscopic redshift l<z<3. In order to have a homogeneous wavelength 
coverage in the spectra, the sample was divided into two sub-groups: galaxies at 1 <z< 1.6 whose spectra cover the rest-frame 
range ~2700-4300A, where the [OII]/13727 emission line can be observed, and galaxies at 1.6<z<3 whose spectra cover the range 
~1100-2800A. In the selection procedure we fully exploit the available spectroscopic information. In particular, on the basis of 
three continuum indices, we are able to identify and exclude from the sample galaxies in which old stellar populations might bring a 
non-negligible contribution to IR luminosity (L /s ) and continuum reddening. Using Spitzer-MIPS and Herschel-PACS data we derive 
L[ R for two-thirds of our sample. The Lut/Lyy ratio is used as a probe of effective attenuation ( Ai RX ) to search for correlations with 
continuum and spectroscopic features in order to derive empirical calibrations to correct UV and [OII]/i3727 luminosities for dust 
extinction. 

Results. Through the analyses of the correlations between different dust attenuation probes, a set of relations is provided that allows 
the recovery of the total unattenuated SFR for star-forming galaxies at l<z<3 using UV and [OII],13727 luminosities. 

The relation between A IRX and UV continuum slope (J3) was tested for our sample and found to be broadly consistent with the literature 
results at the same redshift, though with a larger dispersion with respect to UV-selected samples. 

We find a correlation between the rest-frame equivalent width of the [Oil],(3727 line and / 3 , which is the main result of this work. 
We therefore propose the rest-frame equivalent width of the [OII],i3727 line as a dust attenuation probe and calibrate it through A IRX , 
though the assumption of a reddening curve is still needed to derive the actual attenuation towards the [OII],l3727 line (A[ 0/ ,]). We 
tested the issue of differential attenuation towards stellar continuum and nebular emission: our results are in line with the traditional 
prescription of extra attenuation towards nebular lines. 

Finally, we use our set of cross-calibrated SFR estimates to look at the relation between SFR and stellar mass. The galaxies in our 
sample show a close linear relation (cr=0.3 dex) at all redshifts with a slope ~0.7 - 0.8, which confirms several previous results. 

Key words, galaxies: star formation - galaxies: high-redshift - dust, extinction - infrared: galaxies - galaxies: evolution - cosmology: 
observations 


1. Introduction 


In the past decades, great effort has been devoted to the study of 
the Universe in the redshift range 1 <z<3. Several studies have 
shown that this is the epoch when a substantial fraction of galaxy 
mass assembly took place, and when there is a peak in the evo- 
lutio n of the star formation rate (SFR) density through cosmic 


time dLil lv et alJl996UMadau et al.l l996t Di ckinson et al.l2003l: 


iHopkins & Beacoml 120061 : iDaddi et alJ 12007 ). In this epoch a 
critical transformation phase is believed to have occurred, re¬ 
vealed by the observed changes in the colour-mass plane where 
a significant fraction of galaxies moves from the blue cloud of 
active star formation to the red sequence inhabited by spheroidal 
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galaxies with weak or suppressed sta r formation at z<l (e.g. 
ICassata et al .1120081: ICimatti et al.ll2013l) . 

A key element in the study of galaxy evolution is obviously 
the availability of large samples of galaxies. Spectroscopic red¬ 
shift surveys play a crucial role as they provide samples with 
confirmed redshifts. Though photometric redshift surveys have 
now reached a high level of accuracy, the number of catastrophic 
failures, even if low at a few percent (Ilbert et al. 2013), could 
still produce large unknowns. Moreover, spectroscopy provides 
information that is not accessible by broad-band photometry, e.g. 
emission line fluxes and absorption features that are fundamental 
to gaining insight on gas properties and therefore dust extinction, 
star formation, and feedback mechanisms. 

Astronomical observation plans are actually going in 
the direction of collecting spectroscopic information for 
increasingly large samples of galaxies, as confirmed by 
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the number of spectroscopic surveys that have been car¬ 
ried out in recent years or are planned for the n ear fu¬ 


ture (e.g. K20 jMignohetalJ 2005|), VVDS dLe Fevre et al 


20051) . Deep2 (IWillmer et alJ l2006h. zCOSMOS dLilly et al] 
2007), ESQ-GOODS (IVanzella etak 120081: iPopesso et al]|2009t 


Balestra et al] l20lob . PRIMUS "(iCoil et al.l 1201 ll) . GMASS 


dKurk et ah 20131) . DEIMOS/Keck (Kartaltepe et al., in prep.; 


Capak et al., in prep.), VUDS (iFe Fevre et al.ll2015b. MO S DEF 
dKriek et all l20l5l). VANDEL^l BigBOSS (ISchlegel et al 


20111) . Euclid dLaureiis etaD 1201 lh . WFIRST dSnergel et al 


,201 3)). The large number of high-redshift galaxy spectra that 


is rapidly becoming available inspires the need to test existing 
analysis methods and calibrations, and to create new ones, to 
correctly interpret the wealth of spectroscopic information in the 
frame of understanding the main processes that regulate galaxy 
evolution. 

The rate at which a galaxy produces stars is a critical ingredi¬ 
ent to the investigation of galaxy evolution. Though other prop¬ 
erties like stellar mass can be quite robustly estimated from SED 
fitting to photometric data, there is still room for improvement 
in our ability to estimate the star formation rate (SFR) at high 
redshift because of the uncertainties in the assumptions made 
in tra nslating integrated fluxes and colours to an estim ate of the 
SFR dMadau & Dickinsx)nll20l4tlPannella et alJl2QI5l) . 

The rest-frame UV continuum light emitted by young mas¬ 
sive stars and the Hu optical emission line are _very good 
and primary tracers of star formation (e.g. IKennicutt! Il998t 
iKennicutt & Evansl2012l) . At high redshift, where Ha is less eas¬ 
ily accessible, the [OII]/l3727 emission line can be a fair alterna¬ 
tive. However, these primary and secondary tracers are known to 
be affected by the presence of dust that absorbs the flux emitted 
at UV and optical wavelengths, and re-emits it at longer wave¬ 
lengths in the infrared (IR) regime. The use of dust-unbiased 
tracers like mid- and far-IR (FIR) or radio continuum would be 
preferable, but it is still limited at high redshift due to sensitivity 
limits. Giant strides have been made in this respect with the ad- 


e.g. Futz et al. 2011 

Rodighiero et alJl2010allElbaz et al. 201 ll 

iNordon et alJ 2013: 

Oteo et al. 2014), but the low-SFR regimes 

at high redshift can 

still only be accessed by means of stacking 


(iRodighiero et al.ll2014l ; lPannella et aljl2015l) . In order to study 
star-forming galaxies in a large range of SFR, UV and line emis¬ 
sion luminosities are still the best choice for cosmologically rel¬ 
evant galaxy samples, and the study of reliable corrections for 
dust attenuation continues to be a crucial topic. 

Dust attenuation affecting emission line luminosities can be 
derived by measuring the Baimer decrem ent in galaxy spectra 
dBrinchmann et al Jl2004tlGarn et al.ll201(]h . but this information 
is still rarely available at high redshift. Indirect ways to cor¬ 
rect for dust extinction are often employed, such as the com¬ 
parison with other estimates of SFR like the value derived from 
SED fitting to broad-band photometry (iForster Schreiber et alJ 
20091) u nder the assumption of some attenuation curves (e.g. 


Calzetti et alJl2000l) . 


The dust correction of UV light commonly relies on 


uum 

and dust attenuation (Meurer et al. 19991 

1 Calzetti et al. 

2000 

; DaddietaL 20M Overzier et al. 20ljF 

Takeuchi et al. 

2012 

), though the general validity of such correlation, espe- 


cially at high redshift, is still under debate (ICalzettil 2001 


Boissier et al.l2007l:ISeibert et al.l2005l Reddv et alJ2010L 2012 . 

2015t Overzier et al. 2011 : Takeuchi et al. 2012t Buat et al 


2012; HeinHetak 2013|_ Nordonetak 2013; ICastellano et al.l 
2014tlQteo et al.ll2013112014 lHathi et al.ll2015l) . 


Though the data coming from recent IR surveys are not suffi¬ 
cient to map the population of SFGs down to low-SFR regimes, 
they can be effectively used to test the existing calibrations of 
dust extinction correction recipes, or to derive new ones. This 
can be made under the assumption that the IR emission depends 
only on the absorption of the flux emitted by the young stars re¬ 
sponsible of UV light. However, there might be a non-negligible 
contribution of old stars to the dust heating that grows with de¬ 
creasing sSFR (Id a Cunha et alJ 20081: IKennicutt & Evansll2012l : 


I Arnouts et al.l 1201 3b lUtomo et alJl2014t) . leading to an overcor- 
rection of UV light when IR luminosity (L/«) is used as dust 
correction probe. A way to overcome this problem is to use 
spectroscopic samples, where continuum indices may be com¬ 
puted that indicate the possible presence of old stellar popul a- 
tions dBruzual All 1 98.31 iDaddi et aPl2005UCimatti et al.ll2008l) . 


In this paper, different indicators of star formation for a sam¬ 
ple of 1 < z < 3 galaxies will be analysed. We concentrate on a 
well-controlled spectroscopic sample, rich of ancillary panchro¬ 
matic photometric data, including IR from Spitzer and Herschel, 
in order to focus on the information that can be obtained from 
spectra. We will use the IR-derived SFR as a benchmark against 
which compare other estimates at shorter wavelengths, namely 
UV and [OII]/13727 luminosities, and to derive new (or update 
existing) calibrations of dust extinction correction recipes. This 
way we are able to build a set of self-consistent recipes to derive 
the total un-extincted SFR, which are collected in a Table at the 
end of the paper. 

Throughout this paper, we adopt Hq-1Q km/s/Mpc, Q m =0.3, 
O,\=0.7, give magnit udes in AB photometric system, and as¬ 
sume a lKrounal(l2Q0ll) initial mass function. All linear fits, unless 
differently stated in the text, will be ordinary least squares (OLS ) 
regressions of Y on X (Y|X) feevingtonll969l : llsobeet al.ll990b . 
Finally, we will indicate with the suffix ”0” all un-extincted 
quantities, i.e. quantities corrected for dust extinction. 


2. The multi-wavelength dataset 

The main ingredient of the present study is a well-controlled 
spectroscopic sample at 1 < z < 3 selected from the pho¬ 
tometric catalogue of the Galaxy Mass Assembly ultra-deep 
Spectroscopic Survey (GMASS). 


2.1. The GMASS survey 

The GMASS0 survey dKurk et al.ll2013l) is an ESO VET large 
program project based on data acquired using the FOcal Reducer 
and low dispersion Spectrograph (FORS2). The project’s main 
science driver is to use ultra-deep optical spectroscopy to mea¬ 
sure the physical properties of galaxies at redshifts 1.5<z<3. The 
GMASS photometric catalogue is a pure magnitude limited cat¬ 
alogue from the GOODS-South public image in the IRAC band 
at 4.5 pm: »?4 s<23.0 (AB system). In the chosen redshift range 
this selection is most sensitive to stellar mass. In particular, the 
limiting mass sensitivities are log(M/M Q ) ~ 9.8, 10.1, and 10.5 
for z = 1.4, 2, and 3, respectively. The GMASS photometric cat¬ 
alogue gathers information from U band to IRAC 8.0 pm band 
of 1277 objects. 


1 http:/lvandels.inaf.it/ 


2 http://www.mpe.mpg.de/~kurk/gmass/ 
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150 [ 1 1 1 1 | 1 1 1 1 | 1 1 1 1 | 1 1 1 1 | 1 1 1 

Total GMASS sample at z SI (#758) _ 
z spec at z SI (#385) 

_ I SF+mixed galaxies at z spec SI (#327)_ 

SFGs at z sp „ SI (#286) 

jOO - 24,u.m-detected _ 



Redshift 

Fig. 1 . Redshift distribution of total GMASS sample at z > 
1 with either spectroscopic or photometric redshift (grey); 
GMASS sample at z sp ec > 1 with reliable spectroscopic red¬ 
shift (green); preliminary sample of 327 SF+mixed galaxies 
(salmon); final sample of 286 SFGs analysed in the paper 
(black); 24//m-detected SFGs (blue); PACS-detected SFGs (red). 


2.2. Spectroscopic data 

From the GMASS photometric catalogue, 170 objects were 
selected as main targets to be observed (250 objects includ¬ 
ing fillers) during 145h, and a secure spectroscopic redshift 
could be determined for 131 of them. The spectral resolu¬ 
tion of GMASS spectra is R — A/AA ~ 600. To extend 
our analysis, we collected public spectra from other spectro¬ 
scopic surveys for the galaxies in the GMASS photometric cat¬ 
alogue with no GMASS spectrum. In particular, we searched 
for counterparts of GMASS galaxies from the following sur¬ 
veys: the ESO-GOODS/FO RS2 v3.0 (IVanzella et alJl2008h and 
ESO- GOODS/VIMOS 2. 0 dPonesso et al .1120091 iRalestra et al.l 
|20K)I). the VVDS vl.O dLe Fevre et al.1120051) . and the K20 
(Mig noli et al.ll200~5l) . We refer the reader to the cited papers for 
more information about the various surveys. 

A total of 385 galaxies with a reliable spectroscopic redshift 
Z > 1, coming from any one of the cited spectroscopic surveys, 
was selected (see Fig.Q}, 


2.3. Mid- and far-IR photometry 

The GOODS-South field has been observed with the 
Photode tecto r Array Camera and Spectrometer (PACS; 
iPoglitsch et al.l ( l2010l) ) on board the Herschel Space 
Observatory!^ as part of t wo pro jects: the PACS Evolutionary 
Probe (PEP ; Lutz et all (1201 lh ) and the GOODS-Herschel 
felbaz et al.l 1201 lh programmes. The publicly released PACS 
catalogue produced using as priors the source positions ex¬ 
pected on the basis of a deep Spitzer-MIPS 24 pm catalogue 


dMagnelli et al.ll201ll) was used in this work. We refer the reader 
to lMagnelli et al.l ~ 20131) for all the details about data reduction 
and the construction of imag es and cata logues. The public 
catalogue presented in iMagnelli et al. ( 2013t) has a 24 pm flux 
cut at 20 pJy. We extended it to fainter fluxes with a 24 pm 
source catalogue from Daddi et al. (in prep.) that was created, 
similarly to the IMagnelli et al. ; i201 it) one, using a PSF fitting 
technique at the 3.6 pm source positions as priors. 

GMASS coordinates of the 385 galaxies at 1 <z S7 , CT <3 were 
matched to the sources listed in the 24 pm source catalogues. 
Counterparts were searched within an angular separation of 1.4”. 
After some tests, this angular separation was chosen as the best 
compromise not to lose real counterparts, while minimizing the 
selection of false pairs. 236 24 pm counterparts were found, 121 
of which have at least one PACS detection. 


3. Selection of the star-forming galaxies (SFGs) 

Since the aim of this work is to study the star formation in star¬ 
forming galaxies (SFGs) at 1 <z<3, the catalogue assembled of 
385 galaxies at z jpec > 1 had to be cleaned of quiescent objects 
and active galactic nuclei (AGNs). 

On the basis of the spectroscopic features, 33 quiescent 
galaxies were identified and excluded from the sample. Their av¬ 
erage spectrum (Fig.Q shows the characteristic features of the 
rest-frame UV spectrum of old and passive stellar populations: 
the very red continuum, characterized by prominent breaks (at 
2640A, 2900A, 4 OOOA) and rich of metal absorptions ( e.g. Call 
H&K doublet) (ICimatti et al .12008! iMignoli et alj2005l) . and the 
weakness or complete absence of [OII]/l3727 emission line. We 
will come back on the shape of the continuum of old stellar pop¬ 
ulations later in this section. 

A combination of spectroscopic features and X-ray informa¬ 
tion was instead used to identify AGN hosts. In particular, we 
excluded galaxies showing typical broad and/or narrow emission 
lines (e.g. CIV/11549A, MgILt2800). We also excluded galax- 
ies w ith an X-ray detection from the Chandra 4Ms catalogue 
JXue et al.ll201 ll) . and absorption-corrected rest-frame 0.5-8 kev 
luminosity higher than L x = 3x 10 42 erg/s, following the criteria 
for AGN identification given bv IXue et al] (1201 ll) . This selection 
made us discard 25 more galaxies from the sample. The objects 
with an X-ray detection, but L x < 3xl0 42 erg/s were instead kept 
in the sample. 

The remaining, at this stage, 327 galaxies showed spectro¬ 
scopic features typical of SFGs and lack of AGN evidence from 
the X-rays. UV spectroscopic features identifying a SFG are the 
presence of strong nebular emission lines, such as [OII]/l3727, 
and/or a blue UV continuum. Not all the features can be de¬ 
tected in all the spectra because, given the wide range of redshifts 
covered by our sample, the galaxy spectra sample different rest- 
frame wavelength ranges. In particular, we can roughly divide 
the sample into two sub-groups: galaxies at 1 <z< 1.6 that cover 
the range ~ 2700-4300A ([Oil] sample), where the [OII]/13727 
emission line can be observed; and galaxies at 1,6<z<3 that 
cover the range ~ 1100 - 2800A (UVsample), in whose spectra 
strong inter-stellar medium (ISM) absorption lines can be de¬ 
tected (see Fig.Q. 

3.1. Spotting the presence of old stellar populations in SFGs 


3 Herschel is an ESA space observatory with science instruments pro¬ 
vided by European-led Principal In vestigator co nsortia and with impor¬ 
tant participation from NASA dPilbratt et all20ld ). 


In this work we decided to use L/r as benchmark SFR indi¬ 
cator, exploiting the recent HERSCHEL data in combination 
with deep spectroscopy. L/r is the bolometric luminosity com- 
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Fig. 2. From top to bottom: a) average spectrum of 33 quiescent galaxies in the GMASS sample; b) average spectrum of mixed 
galaxies in the [Oil] sample, at 1 <z< 1.6; c) average spectrum of not -mixed SFGs in the [Oil] sample, at 1 <z< 1.6. Mixed galaxies 
have high values of the three continuum indices indicated by grey shaded regions: MgUV (2600 - 2900A) (iDaddi et alJl2005b : 
C(29 - 33) (2900 - 3300A) (ICimatti et al.ll2008l) : D4000 break (3750 - 4250A) (iBruzual A.II1983I) . See the text and Fig.gjfor more 
details. 


ing from dust emission and can be directly converted into the rate 
of ob scured star formation (iKennicutdl 19981 iKennicutt & Evansl 
120121) . The validity of such conversion relies on the assump¬ 
tion that only young star populations (the ones responsible for 
the UV luminosity emission) heat the dust that than re-emits in 
the IR. However, there might be a non-negligible contribution of 
old stars to the dust heating that grows with decreasing sSFR 
( da_CurthaetaL[ 2008fc_ Kennicutt & Evansl l20l2t lArnouts et al.1 
120 1 3b Utomo et al.l 20141) . leading to an overestimate of the true 
SFR when using as sole estimator. 


Thanks to the spectroscopic information we are able to over¬ 
come this issue. In the wavelength range covered by our spec¬ 
tra three continuum indices may be computed that indicate the 
possible presence of old stellar populations. These indices are: 
MgUV, which is a feature produced by the combination of the 
strongest breaks and absorptions at 2600 - 2900A (iDaddi et al.l 
120051) : C(29 - 33), which is a colour index of the UV continuum 
defined at 2900 - 3300A (ICimatti et al.ll2008h : and the D4000 
break dBruzual A.lll983h . We refer the reader to the cited papers 
for the definitions of the indices. The range over which each pa¬ 
rameter is defined is indicated in Fig. [2] 


Typically, early-type gala xies have MgUV > 1.2, 

C(29 - 33) > 0.6, £>4000 > 1.6 (IDaddi et al.ll2005UCimatti et al 


2008 


[Mignoli et al.l 12005b iKauffmann et al.ll2003t lHathi et al 


2009)- We measured the indices on the spectra of the 327 


galaxies resulting from the previous selection in order to search 
for evidence of the possible presence of old stellar populations 
along the young ones. 


[Oil] sample 

Given the different rest-frame wavelength coverage of the 
spectra, not all indices are available for all galaxies. The spectra 
of the galaxies in the [Oil] sample always have at least one 
index (often two) in their range. We find that 22% of the [Oil] 
sample have a continuum suggesting the presence of an old 
stellar population that might contribute to dust heating and to the 
reddening of the UV continuum. The parameters distributions 
are shown in Fig. [4] and the mean value of each index measured 
on the 33 quiescent galaxies is also indicated as reference. To 
make a distinction, we will call the objects whose indices are 
above the threshold mixed galaxies (because they show both the 
[OII]/13727 line and a red continuum), while the objects whose 
indices are below the thresholds will be referred to simply as 
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SFGs in the rest of the paper. The relative redshift distributions 
of the two classes are shown in Fig. [3 The average spectrum 
of mixed galaxies, presented in Fig. [3 though showing the 
[Oil] /13727 emission line, has a continuum shape more similar 
to that of quiescent galaxies, whose average spectrum is shown 
as reference. Only SFGs (i.e. not -mixed) are used in our analysis. 

UV sample 

In the UV sample only the MgUV index is available in the spec¬ 
trum of ~ 45% of the sample. For the remaining ~ 55% of the 
UV sample no index is available. When in range, the MgUV in¬ 
dex is always < 1.2, within the errors. Moreover, all the galaxies 
in the UV sample (with and without measurable MgUV index) 
have the bluest colour in the (NUV-r) vs. (r-K) colour-_colour 
plot, implying that they have the highest sSFR (lArnouts et al.l 
120131) . This means that in these galaxies a possible contribution 
of old stellarjpopulations to dust heating should be negligible 
(see also lUtomo et al.l2014l) . Therefore we classify all the galax¬ 
ies in the UV sample as not -mixed SFGs, and include them all in 
the analysis. 

3.2. Summary of the final sample selection 

From the magnitude limited GMASS photometric catalogue we 
preliminary selected 327 galaxies with a secure spectroscopic 
redshift 1 <z<3, typical spectroscopic features associated to star 
formation (presence of [OII]/I3727 emission line and/or a blue 
continuum), and lack of evidence of AGN presence from the X- 
rays. We summarize the main properties of the parent spectro- 
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Fig. 4. Continuum indices distribution over the [Oil] sample: 
MgUV, C(29-33), D4000. In each plot the blue vertical line rep¬ 
resents the threshold that separates SFGs from mixed galaxies 
(see the text for more details). The red lines mark the average 
value of each index measured on the GMASS quiescent galaxies 
at Zspec 2 1, for reference. 



Redshift (z ) 

\ spec' 


Fig. 3. Redsfhit distribution. Big plot: the grey histogram repre¬ 
sents the preliminary selection of 327 SF4 -mixed galaxies; the 
blue histogram represents the [Oil] sample (i.e. spectra cover¬ 
ing the wavelength range ~ 2700 - 4300A); the red histogram 
represents the UV sample (i.e. spectra covering the wavelength 
range ~ 1100 - 2800A). Little inset: zoom of the [Oil] sample: 
empty blue histogram represents the total [Oil] sample, while 
light blue shaded histogram represents only not -mixed SFGs in 
the [Oil] sample. 


scopic sample in Table Q] and give additional information on the 
preliminary sample of 327 galaxies in Tablc[2] A deeper spectro¬ 
scopic analysis of the preliminary selection of galaxies revealed 
the possible presence of old stellar populations in ~ 12% of 
them. These 41 mixed galaxies will be excluded from the analy¬ 
sis because of the possible contribution of old stars to dust heat¬ 
ing and continuum reddening. The 286 SFGs will be used to 
study the relations between different estimators of SFR. In the 
rest of the paper, when talking about SFGs we will refer to the 
not -mixed ones. In Table [3] we report some additional spectro¬ 
scopic information about the 286 SFGs. 

On the basis of the wavelength coverage of the spectra, the 
sample can be roughly divided into two sub-groups: galaxies at 
1 <z< 1.6 whose spectra cover the range ~ 2700 - 4300A, and 
galaxies at 1.6<z<3 whose spectra cover the range ~ 1100 - 
2800A. We will refer to these groups, respectively, as the [Oil] 
sample and the UV sample. In Fig. [5] the average spectra of the 
galaxies of the two sub-samples are shown. 

3.3. Comparison to the parent sample 

It is important to establish which part of the galaxy population 
at 1 <z<3 is represented by our sample. The main bias of our 
sample comes from two elements: the requirement of a spec¬ 
troscopic redshift, and the different detection limits of the var¬ 
ious sets of photometric data, in particular the IR ones from 
Spitzer and Herschel. To place our galaxies into context, we 
compared their distributions of mass and colours to those of the 
full GMASS sample in the same redshift range (Fig. [6]). Here we 
note that GMASS is basically a mass-selected sample, whose 
limiting mass sensitivity over the entire redshift range of our 
choice (l<z<3) is log(M/M Q )~ 10.5 dKurk et al.112013l) . 
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Table 1. Summary table of the general properties (spectral classes. X-ray and IR detections) of the GMASS spectroscopic sample at 

Zspec — 1 ■ 


GMASS galaxies with spectroscopic data and reliable spectroscopic redshift ( z spe c > 1) 

[Oil] sample 219 objs. 

quiescent 

SFG 

AGN 

L z >3xl0 42 erg/s 

L z <3xl0 42 erg/s 

NO L x 

24 p m 

PACS" 


quiescent 

21 

0 

0 

0 

4 

17 

6 

4 

Spectral 

SFG 6 


192 

0 

4 

15 

173 

130 

72 

class 

AGN 



6 

2 

0 

4 

3 

1 


L x >3xUy*~erg/s 




6 

0 

0 

6 

3 

X-ray 

L x <3xl0 42 erg/s 





19 

0 

15 

11 


NO L x 






194 

118 

63 

IR 

24 p m 







139 

77 

data 

PACS" 








77 

UV sample 166 objs. 

quiescent 

SFG 

AGN 

L x >3xl0 42 erg/s 

L*<3xl0 42 erg/s 

NOL* 

24 p m 

PACS" 


quiescent 

12 

0 

0 

1 

1 

10 

4 

3 

Spectral 

SFG 6 


148 

0 

9 

3 

136 

88 

36 

class 

AGN 



6 

5 

0 

1 

5 

5 


Lx>3xl0 42 erg/s 




15 

0 

0 

12 

10 

X-ray 

L x <3xl0 42 erg/s 





4 

0 

2 

1 


NO L x 






147 

83 

33 

IR 

24 p m 







97 

44 

data 

PACS" 








44 


Notes. (a> Detection in at least one of the three PACS bands . (b) Preliminary selection of SFGs. It includes both not -mixed SFGs and mixed galaxies. 
See the text for more details (Sec. 13.11 . 

Table 2. Summary table of the general properties (spectral classification and continuum indices. X-ray, and IR detections) of SFGs 
and mixed galaxies (Tot. 327 galaxies). 



z 

Spec. 

class 

Lx 

(erg/s) 

Tot. 

num. 

MgUv" 

C(29-33) fl 

D4000 n 

no IR 

MIPS 
24 p m 

PACS 

TOp m 100/r m 160 p m 

UV sample 

1.6-3.0 

pure SFG 
pure SFG 

<3xl0 42 

NO 

3 

136 

1.07 

- 

- 

1 

58 

2 

78 

0 0 1 

4 25 19 

[ Oil] sample 

1.0-1.6 

pure SFG 
pure SFG 
mixed 
mixed 

<3xl0 42 

NO 

<3xl0 42 

NO 

9 

138 

6 

35 

1.06 

1.49 

0.31 

0.74 

1.30 

1.80 

0 

50 

2 

10 

9 

88 

4 

25 

2 6 5 

11 42 37 

2 4 4 

2 12 9 


Notes. Median values over the sub-sample. The distribution of each parameter in the [Oil] sample is shown in Fig. [4] 


Compared to the reference sample, the SFGs sample is 
probing intermediate-mass objects, the median mass being 
M*=10 9 ' 8 M o , therefore our sample is not mass-complete. The 
sample is biased towards bluer and less massive objects, as 
can be seen from the colour-mass diagram (Fig. [6] c)), with 
the reddest and more massive tail populated mainly by the IR- 
detected galaxies. Objects in the red sequence may be either qui¬ 
escent/passive galaxies or dusty SFGs. The degeneracy can be 
broken by using colour-colour diagrams (Fig. [6] a), b)). The po¬ 
sition of our selection of SFGs in the (NUV-r)vs.(r-K) and the 
(U-V)vs.(V-J) diagrams indicates that we have selected galax¬ 
ies with the highest s SFR, with respect to the parent sample 
dWilliams et al.li2()09t lArnouts et al.lI2013h . of which the IR- 
detected are the dustier ones. 

The different detection limits of the used datasets translate 
into lower limits in the SFRs that can be probed by a specific 
indicator. In particular, UV- and optical-based SFR estimators 
are the deepest ones and allow values down to a few M Q yr~ l 
to be reached at z~ 1. Mid-IR data have slightly higher SFR 
thresholds, while far-IR data are available only for the most star- 
forming galaxies: ~ 100A7 n yr 1 at z~2 an d ~ 10Af o vr 1 at z~l 
dElbaz et al.ll201 it iRodighiero et al.ll201 ll) . 


In the regime of star formation probed by mid- and far-IR 
data we will derive new relations to correct UV- and optical- 
based estimators for dust extinction. Then, we will test the ro¬ 
bustness of extending such calibrations in lower star formation 
regimes by comparing one another the newly calibrated estima¬ 
tors for galaxies without IR data. 

Though our sample is inevitably not complete. Fig. [6] shows 
that it is fairly representative of a specific galaxy population 
whose properties may be summarized as follows: 

- redshift between z ~ 1 and z ~ 3; 

- intermediate stellar mass (approximately in the range 
10 9 - 2 <M/M o <10 10 - 2 ); 

- blue rest-frame colours: (NUV - r )<3 & (r - X)<1.4 or 
(U- V)<1.2&(V-/)<1.5. 

All the relations derived in the paper should be used within 
these boundaries. We also note that the use of our relations is 
not recommended for galaxies whose spectral continuum indices 
show strong evidence of old stellar populations {MgUV> 1.2, 
C(29 - 33)>0.6, D4000>1.6) because they could cause an over¬ 
correction of dust extinction of either UV or [OII]/l3727 lumi- 
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Fig. 5. Top: average spectrum of SFGs in the UV sample, at 1.6<z<3; bottom: average spectrum of SFGs in the [Oil] sample, at 
l<z<1.6. 


nosity and consequently an over-estimate of the total SFR of the 
galaxy. 


Table 3. Number of spectra collected from the different surveys, 
for the final selection of 286 SFGs. The observed spectral reso¬ 
lution (R = A/AA) for each survey is also indicated. 


R 

GMASS 

600 

ESO F oj(5 

660 

ESOv/mos 

180,580 

VVDS 

230 

K20 

260,380,660 

N. 

131 

78 

38 

9 

30 


one Herschel band), we performed an SED fitting procedure, 
based on x 1 minimi zatio n, using the MAGPHYS code in its de¬ 
fault configuration (Ida Cunha et all 120081 ) and all the available 
photometric information, from U band to PACS data. The in¬ 
clusion of the whole SED in the fitting procedure allowed us 
to fully exploit the photometric information. Following an es¬ 
tablished procedure (iRodighiero et ah 201 Oalbt iGruppioni et alJ 


120101 : ICava et all 120101: |Popesso et al. 2012 ). we anchored the 


spectral fit to the FIR Spitzer-Herschel data-points by increasing 
the photometric errors of the bands from U to 8/rm (up to 10%O 
For each source, the 50th percentile of the L/« distribution was 
taken as the best-value; the quoted uncertainty is the 68th per¬ 
centile range (~ 1 cr). 


4. The benchmark: SFR based on IR luminosity 

The total infrared luminosity (L /K ) is defined as the integrated 
luminosity between 8-1000 pm. To derive Ljr for the sources 
with at least two IR photometric points (i.e. 24 pm plus at least 


4 For 8 PACS-detected galaxies we were unable to correctly repro¬ 
duce the observed IR SED. This could indicate the presence of an ob¬ 
scured AGN, whose contribution is unaccounted by our chosen IR SED 
models. To be conservative, we decided to exclude these galaxies from 
all further analysis. 
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8 9 10 11 12 

log(Mass) [M 0 ] 

Fig. 6. From top-left, clockwis e: a) rest-frame (NUV-r)vs.(r-K) colour diagram dArnouts et al. l l2013t) : b) rest-frame (U-V)vs.(V-J) 
colour diagram ([Williams et al.lt2 009h : c) rest-frame colour (U-B) vs. stellar mass diagram. The red dotted line is the red sequence 
at z ~ 2 defined in lCassata et al.l (120081) . In all plots, small light grey dots represent the GMASS parent sample at z> 1, while large 
points represent the final sample of 286 SFGs analysed in the paper. In particular: dark grey points indicate galaxies without IR data; 
salmon points indicate galaxies with a 24 pm detection; red points indicate galaxies with at least one PACS detection. The distribution 
of each quantity is shown by light grey histograms for the GMASS parent sample (#758), while dark grey histograms represent the 
distribution over the SFGs sample; salmon histograms represent galaxies with a 24 /jm detection; red histograms represent galaxies 
with at least one PACS detection. 


To assess the robustness of the resulting quantity, Lm was 
also computed by fitting the mid- to far-IR ph otome try alone to a 
compilation of IR SEP template libraries ( Chary_&Elbaj |200 It 
iDale & Heloull2()0l iLapache et al.ll200fl iRieke et al. 1120091) 


ing the LePhare code ( Arnouts et al.l Il999b lllbert et al.l 20061) . 

The two estimates were found to be highly consistent (s ee Fig . 
0. This test was in line with the results by iBerta et al.l (|20 1 3l) 
who demonstrated that the different methods commonly used in 


the literature to recover L/* all give consistent estimates when 
Herschel photometry is implemented in SED fitting. 

To improve the statistics of the sample, especially at low 
luminosities, L/« was computed also for those galaxies with a 
detection at 24 pm but no PACS data. At high redshift (z~2), 
locally calibrated IR templates have difficulties in recovering 
bolometric Ljr of highly luminous galaxies fr om the 24 pm 
flux alone (lElbaz et alll201 lblNordon et al.ll2010h . Several meth- 
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Fig. 7. L/it computed using UV to far-IR photometry and the 
MAGPHYS code vs. L IR computed using Mid- to far-IR pho¬ 
tometry and the LePhare code. Points are colour-coded with re¬ 
spect to their redshift. Correlation coefficient r xy (logarithmic 
scale) is also indicated. 


ods have been proposed to re-calibrate the 24 nm flux de nsity 
dMurohv et al.ll201 it lElbaz et akll201 it iNordon et al .1120 1 it) that 
make use of the traditional IR templates whose resulting Lir val¬ 
ues need to be corrected using various prescriptions. The major¬ 
ity of these conversions from 24 /j m to Lj R either produce trends 
and biases with the true Lir (defined as the one derived from 
HERSCHEL data), or include non-linear scaling that strongly 
increase the scatter dNordon et al.ll2013h . 

A different approach is to build brand new SED templates 
taking advantage of the most recent far-IR and sub-mm data. 
This latter method is preferable since it does not need to ap¬ 
ply a posteriori corrections to the resulting Lir v alues . Here we 
adopted the main sequence templates of lMagdis et al.l d2012l) to 
extrapolate Lir from the 24 /.an flux densitiefl We tested this 
choice for the sub-sample of PACS-detected galaxies trough the 
comparison between Lir derived from 24 \im alone and Lj R com¬ 
puted using all available photometric data from U-band to PACS, 
as explained in the previous paragraph. The result of this test is 
shown in Fig. [8] where we also show, for comp arison, the out- 
put of using th e popu lar ICharv & Elba j d200lb templates. The 
ICharv & Elbazl (1200 ll) templates work well up to z~ 1, but at 
higher redshifts they overestimate the true Lir, mainly because 
at such redshifts polycyclic aromatic hydrocarbon (PAH) emis¬ 
sions contribute significantly to the 24 uin band flux a nd are dif¬ 
ficult to mode l (lElbaz et al][201 lHNordon et al.ll2010l) . With the 
iMagdis et al.l (1201 2l) templates the resulting 24 /j.m -derived L ]R 
is much more consistent with our true estimate of total Lir at 
all luminosities, though there is still a slight overestimate of the 
total Lir at z> 1.5, that we attribute to the PAH contribution to 
broad-band flux. We quote the mean scatter from the 1-to-l re- 


5 We did not use the MAGPHYS code to derive L IR for galaxies with 
only a 24 jxm detection (and no PACS data) because the code requires 
IR data at longer wavelengths to reliably estimate L tR . 


lation (0.2 dex) as the error on the estimate of Lir when only the 
24 /jm band is available. 

In the rest of the paper, when naming Lir we will be refer¬ 
ring to the Lir estimated from SED fitting to U-to-PACS pho¬ 
tometry using the MAGPHYS code for galaxies with at least 
two IR photometric points (i.e. 24 /Jin plus at least one Herschel 
band), while we will be referring to the Lir estimated from the 
[Magdis et al. ( 2012 ) templates for galaxies with only a detection 
at 24 fjm but no PACS data. Fig.[9]shows the final Lir distribution 
for the sample of SFGs. 




Fig. 8. Comparison between Lir computed with MAGPHYS us¬ 
ing all the available photometric bands (from U-band to PACS), 
and Lir compu t ed usi ng only the 24 /urn b and and, respectively , 
ICharv & Elbazl d200ll) templates (top) and IMagdis et al.l (12012l) 
templates (bottom). Points are colour-coded with respect to their 
redshift (see Fig.|7]for the legenda). 
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Fig. 9. Distribution of L/r over the sample of SFGs with at least 
a 24 pm detection (black histogram). The red histogram shows 
the L/r distribution of PACS-detected sources only. 


24yU.m-detect.ed SFGs 
7 PACS-detected SFGs 



Thejnfrared luminosities can be converted into SFR using 
the lKennicutt) (1 19981) relatiorQ: 

SFR IR (M Q yr l ) = 2.9 X 10~^ L IR [erg/s), (1) 

In star-forming galaxies SFRir is expected to be a good esti¬ 
mate of the total SFR. However, there might be a fraction of 
UV radiation that escapes dust absorption and that is not incor¬ 
porated into the calibrations. This missing unattenuated com¬ 
ponent varies from essentially zero in dusty starburst galaxies 
to nearly 100% in dust-poor dwa rf galaxies and metal-poo r re- 
gions of more-massive galaxies (IKennicutt & Evansll2012r ). To 
overcome this bias many authors suggest adding another term to 
SFRir to account for the unobscured UV light dPapo vich et al. 
2007t[Rodighiero^t^lj|201T)a H Wuvts et al.ll201 it iNordon et all 

20 13t IKennicutt & Evansll2012l) . This second term is simply the 

SFR derived from UV luminosit\Q not corrected for dust extinc¬ 
tion. We adopted the appropriate relation from lKennicuttl (Il998l) . 

SFR uv (M B yr- 1 ) = 0.9 x lO^L^lSOOA)^^ 1 //^ 1 ), (2) 

and finally define 

S FRjr+i/v-S FRir + S FRuv (3) 

as to tal SFR (obsc ured plus un-obscured) (IWuvts et al .11201 ll 
INordon et al]l2013h . 

6 The original relation assumes a Salpeter IMF. We scaled it by a 
factor of 1.55 to match our chosen IMF (Kroupair2001h . 

7 L v ( 1500A) of galaxies in the UV sample was measured directly 
on rest-frame spectra, corrected for slit-losses. In the [Oil] sample the 
range covered by spectra is shifted towards longer wavelengths: in these 
cases, L v ( 1500A) was measured by computing the absolute magnitude 
in a squared filter X (300A wide, centred at 1500A), using the ob¬ 
served magnitude in the filter Y, which is chosen to be the closest to 
T(X) x (1 + z). 


In Fig. [TO] we plot the total SFR ( S FRir+uv) vs. the IR com¬ 
ponent alone. The contribution of S FRuv to the total SFR is al¬ 
most negligible in ultra-luminous (ULIRGs: L;r/L q >10 12 ) and 
the majority of luminous (LIRGs: 10 u <L/r/Lo<10 12 ) IR galax¬ 
ies, while in normal SFGs ( Ljr/L o <10 u ) its contribution to the 
total SFR grows to about 30%. However, we have tested how ne¬ 
glecting this contribution would alter our results and we found 
that none of the calibrations derived in this paper would signif¬ 
icantly change if S FRjr were assumed as total SFR, instead of 
S FRir+uv- 

In the following sections we shall compare S FRir+uv to the 
other SFR estimators available for our sample, with particular 
emphasis on spectroscopic data. Using the IR data to calibrate 
UV and optical estimators, we can derive a set of relations that 
can be used to consistently compare galaxies in a wide range 
of redshifts (l<z<3) and with different kinds of available data 
(from observed optical photometry and spectra, to IR data). 


5. Dust-extinction correction: SFR from UV 
continuum luminosity 

The UV luminosity emitted by young stars is a direct indicator 
of ongoing star formation. However, it is typically severely ex- 
tincted by the dust surrounding star-forming regions, therefore 
the measured rest-frame UV luminosity must be corrected be¬ 
fore being converted into SFR. Empirical relations are usually 
applied for this purpose, most of them calibrated on local galax¬ 
ies. 




0 0.5 1 1.5 2 2.5 

log(SFR IR ) [Mgyr-i] 


Fig. 10. Main plot: comparison between S FRir+uv and SFRir 
for SFGs with at least one IR detection. Filled points mark 
PACS-detected galaxies. Bottom plot: S FRir+uv over .S' FRir ra¬ 
tio as a function of SFRir. Median ratio in bins of S FRir are 
also shown. Upper plot: black histogram is the distribution of 
S FRir over the entire sample of SFGs; red empty histogram is 
the distribution of PACS-detected sources only. 
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5.1. UV sample: SFGs at 1.6<z<3 


In galaxies dominated by a young stellar population, the 
shape of the UV continuum can be fairly accurately approxi¬ 
mated by a power law F^ccAP, where F A is the observed flux 
(erg s " 1 cnT 2 A -1 ) and P is the continuum slope (ICalzetti et al.1 
1 19941) . iMeurer et al.l (1 19991) found, for a sample of local star- 
burst galaxies, a correlation between the UV spectral slope (J3) 
and the ratio of L/r over UV luminosity (L uv = vL v ( 1500A)), 
not corrected f or extinction, which c an b e translated into atten- 
uation A 1600 A . lOverzier et al.l (1201 lh andlTakeuchi et al.l (l2012h 


independently revised the original IMeurer et al.l (119991) relation 


correcting for the effect of the small IUE aperture used in the 
original calibration, thanks to GALEX data, and using an up¬ 
date estimate for L/r (total IR emission instead of IRAS far- 
IR (FIR) emission in the range 40 - 120 pm). The attenuation 
vs. p relation has since been used for galaxies at various red- 


suits about its general validity dCalzett 

i 2001; Boissier et al. 

2007 

; Seibert et a 

.1 2005| Reddv et al 

201(4 Overziei 

et al. 

2011 

Reddv et al 

20121 Takeuchi et al. 

ZOlJPBuat et al. 

2012; 

Heinis et al. 2013 

Nordonetal. 2013: 

Castellano et al. 

2014; 

Oteo et al. 2013L 2014). Here we tested the relation for our sam- 


pie using L/r derived from our mid- and far-IR photometry, in 
order to obtain an estimate of SFR from dust-corrected UV flux 
consistent with the one from IR+UV defined in the previous sec¬ 
tion. 

The attenuation is defined from the ratio between L/r and un¬ 
corrected Luv (IRX = Lir/Luv). However, the exact definition 
is not unique across the literature (IMeurer e t al.lll999k Buat et all 
120051) . Here we adopt the definition (iNordon et al.l2013 ) 


log(S FR ir+uv ) - log(S FRuy ) + 0.4 Ai RX (4) 

(5) 


a,rx - 2 - 5io8 [Wt~ + 11 


t UV 
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where S FRuv is the uncorrected quantity, as defined in Sec. [4) 
and S FRir + uv is the total SFR. This definiti on is consistent with 
the on es adopted by lOverzier et al.l (1201 ll) and iTakeuchi et alJ 
( 1201 21 ) . 

The slope ft is defined in the wavelength range from ~1250A 
to -2600A dCalzetti et all 19941) . and can be measured either di¬ 
rectly on spectra or by using photometry. In our sample, only 
galaxies in the UV sample have spectra with the right wavelength 
coverage; for the [Oil] sample photometry has to be used. 

The 62 IR-detected galaxies with a measurable spectral slope 
were then used to test the Air X vs. [i relation at 1.6<z<3 and 
the result is shown in Fig. [TT] The continuum slope of each 
spectrum was derived through a linear fit_of the average fluxes 
in the spectral windows defined by ICalzetti et al.l ( 1994 ). in the 
log(F ,\) - log(A) plane. More details about the spectroscopic de¬ 
termination of fl can be found in the Appendix. 

We then derive the best-fit relation by performing a linear fit 
of the form A/rx = Co x p + C\ to the data, 

A/rx = (1.10 ± 0.23) xp spec + (3.33 ±0.24) (6) 

with an r.m.s. dispersion of the residuals ~0.7. 

Our relation [Co=1.10) is broadly consistent with the one by 
iMurnhv et al.l (1201 ll) for a sample of 24/rm-detected galaxies at 
z ~2 (Co =0.75). On the other hand, it is flatter than in the UV- 
selecte d samples of lOverzier et al.l (1201 ll) (local) and lBuat et al.l 
(120121) (z~ 1.5) (respectively, Co- 1.96 and Co=1.70). Also, the 
standard deviation in the latter cases is of the order of ~0.3 dex. 


Fig. 11. A/rx vs. [I for IR-detected galaxies in the UV sample. 
Top: p is measured from the spectra. The black line is a fit to the 
data (Eq.6). Bottom: p is measured from photometry. The red 
line is a fit to the data (Eq.7). The black line is the relation for 

Pspec• 


while our relation is more dispersed. We ascribe it mainly to the 
different sample selections. A low dispersion is expected when 
galaxies are UV-selected local starbursts or high-redshift equiva¬ 
lents, but deviations from a tight relation have been reported both 
for dusty IR-luminous galaxies and for quiescen tly star-forming 
galaxies (Buat et al.ll2012t ITakeuchi et al.lf2()Ij . and references 
therein). Our parent catalogue is mass-selected, therefore our 
sample is expected to include objects with a wide range of star 
formation histories and different levels of starburst activity that 
may be characterized by different attenuation laws and increase 
the scatter in the relation, with respect to a UV selection. In par¬ 
ticular, we found that the objects that deviate the most from our 
relation (i.e. with the highest A/rx at fixed P) have the highest 
Eir. 

The continuum slope can be derived also from photo¬ 
metric data, either by measuring p directly from the best-fit 
SED tem plate to the overall photometric data of the galaxy 
(e.g. lOteo et ahl 120131 . 2014 ). or by fitting the photometric 
points which sample the UV continuum at the redshift of the 
galaxy with a power-law function (e.g. iBouwens et al. 2009 : 
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iBuat et aI1l2012b [Nordon et alJl2013t iPannella et alJl2015h . We 
decided to use the latter method, to avoid model dependen¬ 
cies. Uncertainties were computed by randomizing the observed 
fluxes according to their errors and repeating the fit to derive 
a new slope. Since even with the photometry we cannot probe 
wavelengths as short as ~1200A at z~l, for the photometric de¬ 
termination of the continuum slope we decided to use, at all 
redshifts, a shorter rest-frame wavelength range to derive f3 p hot, 
starting from A rest ~ 1500A. This wavelength coverage is also the 
most popular in the literature, when the slo pe is derived from 


photometry and not from spectra (see e.g. 

Bouwens et al. 20091 

l2012tiNordon et al. 2013tlCastellano et al. 

201 i. 


We computed y 3 p hot in the UV sample for the same galaxies 
for which f3 spec had been already derived and a similar relation 
was found between /3 p i wt and Air X , 


A irx = (1.45 + 0.25) xfiphot + (3.63 + 0.26), (7) 

with an r.m.s. dispersion of the residuals ~0.6. This relation 
gives attenuation values that are consistent with those derived 
from spectroscopy. This result will justify us to use the photo¬ 
metric slope to derive a similar relation for the galaxies in the 
[Oil] sample, where the spectra do not cover the wavelength 
range needed to compute / 3. 

5.2. A little digression about regression methods 

Since our aim is to provide a relation to predict the value of a 
variable from the measurement of another, we adopted an OLS 
(Y|X) fit in the derivation of the A; RX vs. [I relation. The OLS 
(Y|X) fit is the most widely used in the literature and we are 
assuming that ”OLS” and ’’linear regression”, in the absence of 
further specification, means OLS (Y|X) in all the works whose 
calibrations we cite here to compare with our findings. 

A fit that treats symmetrically the two variables is instead to 
be preferred when the goal is to estimate the unde rlying func¬ 
tional relation between the variables (llsobe et alJ 1 1990h . The 
slope of the Air X vs. [5 relation (Co) depends on the attenua¬ 
tion curve. The Co parameter expected from a given attenuation 
curve should be compared to the value derived from a symmet¬ 
ric fit, instead of an OLS (Y|X). We report that an OLS bisector 
fit (llsobe et all 19901) to the data analysed in Fig.QT](top) would 
produce the following relation (see Eq. 6 for the corresponding 
OLS (Y|X) fit) 

Airx = (1.96 ± 0.23) Xy3,pec + (4.22 ± 0.19) [OLSBisector], 

( 8 ) 

while in the case of y 3 p hoi (Fig.QT]bottom) we would have 

Airx - (2.16 ± 0.27) x/3 p i lo , + (4.36 + 0.27) [OLS Bisector], 

(9) 

whose corresponding OLS (Y|X) fit is given in Eq. 7. For com¬ 
parison, the Co value expected from the Calzetti law is Co=2.3, 
which gives a relation fairly consistent with our findings. 

5.3. Extinction from observed colour 

Another popular empirical recipe to derive a dust extinction cor¬ 
rection to be applied to UV fluxes, for galaxies at z>1.4, comes 
from the tight correlation that lDaddi et al.l (2004 ) found between 
the observed ( B—z ) colour and E(B-V) of BzK galaxies. This re¬ 
lation is valid only for galaxies in the redshift range 1.4<z<2.5, 
where the ( B - z) colour is basically another way of express¬ 
ing the UV continuum slope. In the iDaddi et al.l (2004) paper. 


the E(B-V) is derived from SED fi tting to UV-to-NIR (near- 
IR) multi-band photometry (assuming lBruzual & Chariot] (12 0031) 
models, Salpeter IMF, constant star formation, and a Calzetti 
law). 



Fig. 12. Airx vs. the (B- z) co lour. Black line is a fit to the data. 
The magenta line is the IDaddi etliil (120041) relation, plotted as 
reference. 

In this work we use instead the available IR data to calibrate 
the ( B - z) colour as an extinction estimator. A linear fit to the 
data produces the following relation: 

A,r X = (1.64 ± 0.13) x{B - z)ab + (0.77 + 0.13). (10) 

In Fi_g.lT2l we nlot all IR-detected SFGs at 1.4<z<2.5. The orig¬ 
inal IDaddi etakl d2004l) relation is shown as a re fere nccQ. The 
two relations are in general quite consistent, but ours is flatter 
for galaxies with reddest (B-z). 

The attenuation derived from (B-z) colour is consistent with 
the one derived from the continuum slope in the previous sec¬ 
tion, not only for the IR-detected galaxies (that were actually 
used to derive both sets of calibrations), but also for galaxies 
without IR information. 

5.4. [Oil]sample: SFGs at l<z<1.6 

In our sample, only galaxies in the UV sample have spectra with 
the right wavelength coverage to compute the continuum slope. 
However, photometry can also be used for the same task. For 
the galaxies in the [Oil] sample we computed the slope /3 p i lot 
between ~ 1500A and ~2600A, as explained in the previous sec¬ 
tion, and we found the following relation with Air X , 

Airx = (1.03 + 0.26) x/3 phot + (3.54 + 0.25), (11) 

8 The IDaddi et all d2004l) relation is E(B-V) vs. (B-z). In the plot, 
E(B-V) was converted to Airx assuming a Calzetti law, as in the original 
calibration: E(B-V)=A /R X /ki-m , where ki 50 o is the value of the redden¬ 
ing curve at 1500A dCalzetti et alfeOPOl) . 
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Fig. 13. ,4 / A >x vs. Pphot for IR-detected galaxies in the [Oil] sam¬ 
ple. The blue line is a fit to the data (Eq. 11). The red line is the 
Airx vs. ppi w , relation for the UV sample (Eq. 7). 
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pie. However, the difference is not highly significant, therefore 
we see no evolution of the A IRX vs. [1 relation between z~ 1.3 and 
z~ 2.3. Moreover, dust attenuation of the galaxies in both sam¬ 
ples broadly follows the prediction of the Calzetti law (when 
considering the OLS bisector relations; see Sec. E3, though in 
the low-redshift bin the Calzetti law tends to over-predict the at¬ 
tenuation in objects with the reddes_t_slope, with respect to the 
results of our calibrations (see also lPannella et al.Il2015l) . 

The relations derived in this and the previous sections were 
applied to all the galaxies in our sample, choosing the appropri¬ 
ate one for each object to derive A IRX values from the UV con¬ 
tinuum slope and compute the un-extincted SFR ( S FRuy 0 ). The 
final check to assess how good the calibrations are to recover 
the total SFR comes from the comparison between S FRuy 0 and 
S FRir+uv in Fig.[T4lwhere we can see that the two estimates are 
in good agreement. 


6. Dust-extinction correction: SFR from [0II]/13727 
emission line luminosity 

The next SFR estimator to be considered was the forbidden line 
[Oil] Ail21. The luminosities of forbidden lines are not directly 
coupled to the ionizing luminosity, and their calibration as SFR 
indicators is sensitive to density variations in dust reddening, 
che mical abundance, a nd ionization among star-forming galax¬ 
ies dKennicutd 1998: iJansen et al.1 1200 It iKewlev et all 120041 


IMoustakas et al.1 120061) . The excitation of [OII]/13727 is suffi¬ 
ciently well behaved that it can be calibrated empirically as a 
quantitative SFR tracer through Ha, though the [Oil |/Ha ratio in 
individual galaxies can vary considerably and this is the main un¬ 
certainty related to the [Oil]-derived estimate of SFR. The other 
two major sources of uncertainty are the effects of metallicity 
and dust extinction on the conversion of [OII]/!3727 luminosity 
into SFR. 

Various approaches have been proposed in the lit eratu re 
to derive dust-corrected estimates of SFR[oii] 0 . iKewlev et al.l 
( 2004 1. for example, derive a local S FR[oir\ 0 calibration based 
on the intrinsic (e.g. reddening corrected) [011 1/1372 7 lu minos¬ 
ity and the abundance. The local calibration bv IMoustakas et all 
( 2006 1 is instead parametrized in terms of the B-band luminos¬ 
ity to remove the systematic effects of reddening and metallicity. 
Both calibrations use as reference SFR the estimate from Her lu¬ 
minosity, corrected for dust-extinction using the Balmer decre¬ 
ment. 


log(SFR UVo ) [M 0 yr -1 ] 

6.1. The [OIIJA3727 equivalent width vs. p relation 


Fig. 14. Main plot: comparison between S FRjr + uv and S FRuv 0 
corrected using the Ajrx vs. p relations derived in the paper, /j is 
derived from spectroscopy, when possible, otherwise from pho¬ 
tometry (see the text for more details). Upper plot: black empty 
histogram is the distribution of S FRjjy 0 over the entire sample of 
SFGs; black filled histogram is the distribution of S FRuv 0 over 
the sample of 24 p/H-dctectcd sources (both with and without 
PACS data); red empty histogram is the distribution of PACS- 
detected sources only. 


Here we exploit the Spitzer-Herschel data and use SFRjR+uy 
as reference to calibrate a dust-corrected estimate of SFR from 
[OII]/L3727 luminosity for our [Oil] sample of galaxies at 
1.0<z<1.6. 

Integrated line fluxes of [OII]/l3727 were measured and cor¬ 
rected for slit-l osses using the available photometry. We assume 
the iKennicutti (119981) calibration (scaled to a Kroupa IMF) to 
convert [OII]/13727 luminosity into SFR, 

SFR [On]0 (M Q yr- 1 ) = 0.9 x 10 ~* 1 hosiergs' 1 ), (12) 


with an r.m.s. dispersion of the residuals is ~0.8 (Fig. EE 
Comparing Eq. 7 and Eq. 11 we find that in the [Oil] sample the 
slope of the relation (Co) is slightly flatter than in the UV sam- 

9 The OLS bisector fit gives Co=1.49 


where T[o//] 0 is the intrinsic [OII]/l3727 luminosity, corrected 
for dust extinction. We explored the possibility of calibrating 
a relation between attenuation, derived by our IR data, and 
some property of the [OII]/13727 line, in order to obtain a 
self-consistent way of computing the dust corrected S FR[on] a 
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Fig. 15. Top: /3 p i w , vs. log( EW[on]) rest-frame. We use the convention of positive EW for emission lines. The right axis converts 
Pphot to A ipx using Eq. 11. Black points represent IR-detected galaxies, while grey points represent galaxies with no IR detection. 
The blue dashed line is a linear fit to the points. Bottom: average spectra in two bins of rest-frame EW[oii\. Lines of interest are 
labelled. 


from [OII]/l3727 information alone, without the aid of addi¬ 
tional multi-wavelength photometric data. An anti-correlation 
(r xy - - 0.56) was found between the rest-frame equivalent width 
(EW) of the [OII]/13727 line (i.e. the blended doublet) and the 
UV continuum slope (Fig. US, both for galaxies with and with¬ 
out IR data. A linear fit gives the relation 


fipho, = (-1.35 + 0.20) x log(EW res t) + (0.91+0.30) (13) 


with an r.m.s. dispersion of the residuals ~0.6. Since f5 is related 
to dust attenuation (Eq. 11), the relation with EW implies a re¬ 
lation between EW and dust attenuation for the galaxies in our 


sample. Combining Eq. 11 and Eq. 13 we obtairPI 
A irx = (-1.39 + 0.26) X log(EW rest ) + (4.48 + 0.35), (14) 


1(1 We propose as the relevant relation the one between the EW and 
P , instead of the one directly linking the EW to A tRX because /? is the 
only dust attenuation related quantity that we can directly measure in 
all the galaxies in our sample. In the original definition (Eq. 5), A IRX is 
a directly measured quantity, but it is only available for galaxies with IR 
data. On the other hand, A IRX derived from P pho , (Eq. 11) is available for 
the entire sample, but it is a derived quantity. However, we report that 
a direct fit of A IRX derived from P pho , (Eq. 11) and log(EW rest ) would 
give a relation consistent with Eq. 14. The same result would be ob¬ 
tained with a direct fit of A IRX as in the original definition (Eq. 5) and 
log{EW r est), for IR-detected galaxies. 
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Despite all our efforts to clean the sample from galaxies in which 
older stellar populations might contribute to dust heating (i.e. 
Li if) and continuum reddening, we cannot completely exclude 
a mild contribution that might increase the scatter in the /3 p i w , 
vs. log{EW res t ) relation in the low-EW[o//] part of the sample, 
given the presence of Balmer absorption lines (stronger in A- 
type stars) and Call H&K lines in the low-EW[o//j stacked spec¬ 
trum (Fig. IT5l) . On the other hand, we measure a stronger EW 
of the MgII/12800 ISM absorption line (i.e. the blended doublet) 
in the low-EW[o//j stacked spectrum (Fig. 1 1 6[ i: ISM absorption 
lines have been observed to be stronger in SFGs_with redder /3 
and ch arac terized by higher dust extinction (IShaplev et al .120031 
iTalia et alJl2012h . 


Mgll A2800 



Fig. 16. Composite spectra of galaxies in the [Oil] sample in two 
bins of rest-frame EW[oii ]: zoomed-in image of the MgII/12800 
ISM absorption line. 


6.2. Continuum vs. nebular attenuation 


We cannot use directly A/ KX derived in Eq. 14, i.e. the attenua¬ 
tion towards the continuum, to correct the [OIT] 43727 luminos¬ 
ity, first because of its dependency on wavelength, and second 
because of the extra amount of reddening suffered by nebular 
emission, with res pect to the co ntinuum, claimed by many au¬ 
thors starting from lCalzettif (1997!). The first point means that an 
extinction curve must be assumed. It is known that the effects of 
dust are stronger at shorter wavelengths, therefore correcting a 


luminosity measured at 43727A using an attenuation measured 
at 41500A would produce an over-correction. The extra red¬ 
dening towards nebular lines is an issue still strongly debated. 
In fact, though many a u thors have confirmed the original find¬ 
ings by ICalzettil ( 1997^]_ Calzettietak J 20001) also at high red - 
shift (see e.g. Forster Schreiber et al. 20091 : Wuvts et alJfeoi lh . 
other studies claim a ratio bet ween continuum and nebular at¬ 
tenuation much closer to unity (Reddv et al,l201(A iKashino et alJ 


201.1: iPannella et al.ll201-5i:iPuglisi et all2015h. 


To summarize, the dust extinction corrected SFR is 


SFR [Om0 (M Q yr- 1 ) = 0.9xl0- 41 L [o//] xlO 0 ' 4 ^"/^- 1 ), (15) 


where Lyon] is the observed, uncorrected [OII]/l3727 luminosity. 
If we parametrize the extra reddening as 


E{B - V)neb = E{B - V) cont /f (16) 


and assume a reddening curve (kj), then Ayoii] is defined as 

A [on] = A irx x {k h Ji<i 5 oqa) x / 1 • ( 17 ) 

We use the value of the extinction curve at Her as prescribed by 
iKennicutld 19981) because of the manner in which the [OII]/l3727 
luminosities were cali brate d through the [OII]/Hcr ratio. 

In the UV regime a lCalzetti et al.1(120001) law is a fair assump¬ 
tion for the reddening curve, for our sample. In the optical range 
with our data we are unable to constraint the reddening curve 
and we decided to assume the same ICalzetti et al.l d2000l) law 
to ease the comparison with rec e nt works on differ e ntial atten- 
uation (e.g. IKashino et alJ 120131 : 1 Wuvts et~akll2013| : iPrice et ahl 
120141) . though it is still unclear if this is an appropriate assump - 
tion in the context of high-redshift galaxies dReddv et al.ll2015 ). 

We derive dust extinction corrected S ER\ou\ n assuming dif¬ 
ferent values off and compare the result to S FRyuyy 0 (corrected 
using /3 p hot and Eq. 11) to derive an estimate of the best/-factor. 
The best agreement is found with f— 0.50, which is very close to 
the 0.59 value quantified bv ICalzetti! d 1997H 1 *1 . In Fig. [171 where 
we also plot, for comparison, the fits corresponding to other two 
different choices off (0.75 and 1.0). The small dispersion (~0.2) 
is due to the correlation of the corrections on the two axes, since 
we derive A you] from Airx- Though we cannot draw any conclu¬ 
sive statement, our data appear to be more consistent with the 
need of an extra attenuation towards nebular lines in galaxies at 
high redshift, than with studies that claim an / coefficient much 
closer to unity. 



log(SFR [OII)o ) [M 0 yr-‘] 

Fig. 17. S ER/jv,, vs. S FR[oii] 0 corrected for dust extinction us¬ 
ing, respectively, //,/„„ and Eq . 11, and EWyon] and Eq. 14 as¬ 
suming a lCalzetti etal.1 ( 120001 ) law and f- 0 . 50 . The black line 
is the 1 -to-1 relation that also corresponds also to a linear fit to 
the points. Blue and red lines show the fits with other two dif¬ 
ferent assumptions of /. Blue points are IR-detected galaxies, 
while light blue points are galaxies with no IR detection. 


11 The canonical value quoted by the original works dCalzettill 19971 : 
ICalzetti et alJl200(t) is /=0.44, while we quote /=0.59. The difference 
comes from the fact that in the original papers two different extinction 
curves for the nebular and continuum emission were used, while we 
are assuming the same C alzetti law for both emissions llPannella et all 
l2015tlSteidel et all2014l) 
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We stress that the relevant relation that we derive in this 
work is the one linking p p h 0 t and EWyoip- Then, the actual 
attenuation of the [OIT]T3727 line luminosity depends on the 
(/q^/kjjoq^x/" 1 term in Eq. 17, i.e. on the combination of the 
assumptions about the reddening laws respectively for the stel¬ 
lar continuum and nebular emission, and the best / value that 
we derive by forcing agreement between S FRyoii] 0 and our ref¬ 
erence SFR. With this approach, different assumption s ab out the 
reddening law will produce different estimates off (iPuglisi et al.1 
12015lk but the resulting S FR[oii] 0 values will be always, by con¬ 
struction, consistent with one another. 

As an exercise we derived the / value under a different as¬ 
sumption for the nebular attenuation law. We choose a Galactic 
extinction curv e, consistently with the original pr e scripti on by 
ICalzettil (119971) . in particular the ICardelli et al] (119891) one, 
while for the stella r co ntinuum we continued to assume th e 
ICalzetti et al.1 (l2000h law (iReddv et al.l201.4lSteidel et alJ2014h . 
With this choice of attenuation curves the best agreement be¬ 
tween S FR[on] 0 and S FR\u V \ a is found with /=0.37. 

To conclude this section we summarize the steps to derive 
the dust extinction corrected S FR[oii] 0 ■ 

1. Derive the attenuation towards the continuum (A/rx) from 
the rest-frame EWyoii] and its relation with/ 3 p hot (Eq. 13-14). 

2. Convert A/rx into Ayoii] assuming appropriate reddening 
curves respectively for the stellar continuum and nebular 
emission, and an extra-att enuatio n factor (Eq. 1 7). 

In this work we assume a lCalzetti et al.l d2000t) law for both 
the stellar continuum and nebular emiss ion and apply an /- 
factor /= 0.50. If a lCardelli et al.l (Il989t) Galactic extinction 
curve were instead applied to nebular emission, while keep¬ 
ing the Calzetti law for the stellar continuum, the /-factor to 
be used would be /=0.37. 

3. Derive S FRyoip a by correcting the observed |OII]/13727 line 
l umin osity with Ayon] (Eq. 15); in this work we assume a 
iKennicutfl dl998h calibration. 

Finally, we are aware that our relations are taking the cali¬ 
bration of the [OII]/13727 line as a SFR indicator as a closed box 
because with the data in our possess we cannot tackle the issue 
of how the [OII]/Ha' ratio varies across our sample, as a func¬ 
tion of different gas properties. For this reason our method for 
deriving the S FRyon] 0 should be applied carefully to individual 
galaxies. However, our calibration will be useful for computing 
the statistical star formation and extinction properties of large 
high-redshift galaxy samples, when no other information about 
nebular extinction is available. 


7. A quick look at SFR from SED fitting 


We dedicate a brief section to the SFR estimate derived by fitting 
the galaxy UV-to-NIR broad-band photometry to the spectral en¬ 
ergy distribution (SED) of synthetic stellar populations. For this 
exercise we do not consider the FIR data but use only photom¬ 
etry from U band to IRAC 5.8/rm, to see how the results of our 
calibrations compare to a different dust-corrected SFR estimate 
derived without direct information on dust emission. 

A widespread approach is to use models in which the star 
formation history (SFH) is described by an exponentially de¬ 
clining SFR, though in the last years some studies have in¬ 
troduced increasing star formation histories as a better way 
to describe high-redshift star-forming gal axies, namely at z ~2 
dMaraston et al] 120101: Reddv et al.l l2012h . We choose the set 
of models bv iMarastonf i 20051) and used the HyperZ software 


iBolzonella et al.il2000 ) to perform the SED fitting, star forma¬ 
tion histories were parametrized by exponentials (e~ f ^ T ) with e- 
folding timescale r between 100 Myr and 30 Gyr, plus the case 
of constant SFR. We assumed a Kroupa IMF, fixed solar metal- 
licity, a Calzetti law for dust extinction, and a minimum age 
age m in= 0.09 Gyr. We derived from simulations the errors on 
the quantities estimated using SED fitting techniques. We found 
that errors on both S FR$ ed and stellar masses are, on average, 
~ 30% (Bolzonella, private communication). 
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Fig. 18. Comparison between S FRsed and SFR estimated from 
spectroscopic tracers: UV flux, corrected using/?, for galaxies in 
the UV sample (red points), and [OII]/13727 flux, corrected using 
EW[o//], for galaxies in the [Oil] sample (blue points). 


In Fig. [18] (top) we compare the results from SED fitting, 
under the cited assumptions, to the ones previously obtained 
from our calibrations. In particular, for the UV sample we plot 
S FR(j VlJ corrected using f}, while for the Oil sample we plot 
SFR[oii] 0 corrected as explained at the end of Sec. 16.21 (as- 
suming a Calzetti law for both stellar continuum and nebular 
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emission and applying an extra-attenuation f— 0.50). There is a 
general fair agreement between the two estimates at high SFR, 
but in the low-SFR regime the SED estimate is lower than our 
calibrated values, i.e. lower than S FR/r vuv on which our cali¬ 
brations are based. This tail at low SFRsf. d has been reported 
also by other authors (Wuydsetal. j20 1 lb iReddv et al.l l2012t 
lArnouts et al.ll2013t llJtomo et alJl2014h . according to which this 
difference may have multiple causes. Either IR+UV is overes¬ 
timating the true SFR due to a non-negligible contribution of 
old stars to the dust heating or some of the physical assump¬ 
tions made in the SED fitting procedure do not correctly describe 
the analysed galaxies, causing S FR$ed to underestimate the true 
SFR. In our case, we tend to exclude the former explanation be¬ 
cause we specifically excluded galaxies showing evidence of old 
stellar populations in their spectra. To investigate thoroughly if 
and how degenerate SED fitting parameters (i.e. SFH, age, atten¬ 
uation law, metallicity) should be tuned in order to recover the 
true SFR is beyond the scope of this work. Therefore we limited 
ourselves to a quick test and we changed only one of the main 
assumptions, namely the SFH, by forcing all galaxies to a con¬ 
stant SF (and leaving all the other parameters unchanged). We 
note that, in our sample, when r is left as a free parameter only 
for ~ 20% of the galaxies the best^ 2 is obtained with a constant 
star formation history. The results are shown in Fig. [18] bottom 
plot. Assuming a constant SFH the less star-forming galaxies 
are forced to values of their S FR$ed that are in better agree¬ 
ment with our calibrated estimate with respe ct to the resul ts ob - 
tained by leaving r a free parameter (see also IWuvts et al.ll201 ll) . 
though especially in the [OII]sample SFRsed seems to be sys¬ 
tematically slightly higher than our IR-calibrated values. 



9 10 11 


log(Mass) [M 0 ] 

Fig. 19. SFR vs. Mass in three redshift bins, each spanning ~1.2 
Gyrs. For each galaxy, the plotted SFR estimate is S FRjr+uv 
for IR-detected galaxies and SFR[oir\ 0 or SFRuy 0 for galaxies 
with no IR data respectively in the [Oil] sample and in the UV 
sample. In each panel the three lines are the fits (OLS bisector) 
to the data in the three redshift bins. Different colours identify 
each redshift bin. 


8. An application: the SFR vs. Mass relation 

An obvious application of the calibrations derived in the previ¬ 
ous section is the study of the SFR vs. stellar mass relation. 

T he SFR-M^ relat i on is expected to evolve with red¬ 
shift ( Daddietak||20p7|; iNoeskeetak 2007^Wu ytset_aljl2()l it 


[Whitaker et all201^ Heinis et al.1120131 : iBuat et al.ll2014il . there¬ 

fore we divided our sample into three redshift bins, each span¬ 
ning ~1.2 Gyrs (Fig. IT9ll. The SFR for each galaxy is chosen 
as follows: S FR/r, uv for IR-detected galaxies and the appropri¬ 
ate IR-calibrated SFR estimate (S FR[oii] 0 or S FRuy 0 ) for galax¬ 
ies not detected in IR, depending on the redshift. Stellar masses 
were computed using lMarastonl ( 2003 ) models with the assump¬ 
tions outlined in the previous section and constant SF. 

We find that our galaxies lie on a close linear relation 
with cr~0.3 dex and r xy ~ 0.6 at all redshifts. There are weak 
trends both for the slope and the normalization, with the slope 
becoming steeper and the normalization higher with increas¬ 
ing redshift, though the completeness limits and spectroscopic 
requirements of our sample do not allow a robust qualifica¬ 
tion. However, we note that the dispersion that we find in the 
relation is consistent with that often quoted in the literature 
(see e.g. Daddietak 20071 iNoeske et al.l 120071: IWhitaker et al.l 
1201 2t IBuat et al.l 2014 ). Also, the same trend of the normaliza¬ 
tion with redshift has been reported in previous works, while 
there is not a consensus yet on the evolution o f the slope with 
redsh ift (Elbazetak 20Q7yNoeskeetaT |2007|_ WhitakeretalJ 


12012t iHeinis et alJl2013HBuat et alJl2014t lPannella et al.ll2015l) . 

as an effect of different sample selections, completeness limits 
and chosen SFR estimators. In the redshift range 1 <z<3 litera- 
ture quoted values of the slope range from ~0.7 to ~0.9 (see also 
Wuvts et al. 201 It R odighiero et al.ll201~II : iKashino et aljF2013fc 


Rodighiero et al.ll2QI4l) . with which our findings are consistent. 


9. Conclusions and summary 

In this paper we use a sample of galaxies drawn from the 
GMASS survey to study different SFR estimators. Our aim is 
to use IR data to derive empirical calibrations to correct UV and 
[OH],13727 luminosities for dust extinction. To do this, we con¬ 
centrated on a well-controlled spectroscopic sample, rich of an¬ 
cillary data. 

We started by selecting a sample of SFGs with a spectro¬ 
scopic redshift 1 <z<3 from a pure magnitude-limited parent 
sample (m4.5<23.0). In the chosen redshift range this selection 
is most sensitive to stellar mass. In particular, the limiting mass 
sensitivities are log(M/M Q ) ~ 9.8, 10.1, and 10.5 for z = 1.4, 
2, and 3, respectively. Our sample can be divided into two sub¬ 
samples with homogeneous rest-frame wavelength coverage: 

- galaxies at 1 <z< 1.6 that cover the range ~ 2700 - 4300A 
([Oil] sample), where the [OII],13727 emission line can be 
observed; 

- and galaxies at 1.6<z<3 that cover the range ~ 1100-2800A 
(UV sample), in whose spectra strong ISM absorption lines 
can be detected. 

We excluded from the sample all quiescent galaxies and 
AGNs, identified on the basis of spectroscopic features and X- 
ray luminosity. We also used three continuum indices (MgUV, 
C(29-33), D4000) to clean the resulting preliminary selection of 
SFGs from galaxies showing some evidence of the presence of 
old stellar populations that could bias the correct interpretation 
of the dust- and SFR-related observables. The final sample in¬ 
cludes 286 SFGs, of which on third has IR information coming 
from Spitzer-MIPS and Herschel-PACS, and another one third 
has only a Spitzer-MIPS detection. 
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Table 4. Summary table of all the relations derived in the paper. For each relation, observed and derived quantities are indicated, 
as well as slope and intercept of the linear relation, with associated uncertainties. Relations are given in the form: Derived = 
S lope x Observed + Intercept. 


Observed 

Derived 

Slope 

Intercept 

z 

Notes" 

Sec.* 

ft spec 

A/rx 

1.10 ±0.23 

3.33 ± 0.24 

1.6<z<3.0 

spectral coverage: 1200 - 2600 A 

imiAi 

Pphot 

A/rx 

1.45 ±0.25 

3.63 ± 0.26 

1.6<z<3.0 


TsTT 

(B-z) 

A/rx 

1.64 ±0.13 

0.77 ±0.13 

1.4<z<2.5 


1 5JT 

Pphot 

A/rx 

1.03 ±0.26 

3.54 ± 0.25 

1.0<z<1.6 


1541 

log (Ew [om r 

Pphot 

-1.35 ±0.20 

0.91 ±0.30 

1.0<z<1.6 

[0//]/13727 emission line 

IF 

\og(EW [om ) c 

A ^ 

rt-IRX 

-1.39 ±0.26 

4.48 ± 0.35 

1.0<z<1.6 

[0//J/13727 emission line 

16.11 


Notes. Spectroscopic requirements, when needed. (b) Reference to the section in the paper where the relation has been derived. (cJ Rest-frame. 
The convention of positive EW for emission lines is adopted. Irf) A IRX is the attenuation towards the stellar continuum at T=1500A. The attenuation 
for the [OII]/13727 line is defined as A[oir\=Ai RX 'X.(KH a l Kisooa)*/' _ 1 , where kx is the reddening curve. Assuming a Calzetti law for both the stellar 
continuum and nebular emission we estimate /=0.50. Assuming instead a lCardelli et all i j 19891) law for nebular emission and Calzetti law for the 
stellar continuum we estimate /=0.37 ( Sec. 16.2b . 


Though inevitably not complete, especially due to the spec¬ 
troscopic requirement that introduces a bias towards brighter and 
bluer objects, our sample is fairly representative of a specific 
galaxy population whose properties may be summarized as fol¬ 
lows: 


- redshift between z ~ 1 and z ~ 3; 

- intermediate stellar mass (approximately in the range 
10 9 ' 1 2 <M/M o <10 10 ' 2 ); 

- blue rest-frame colours: (NUV - r )<3 & (r - K)<\A or 
(t/-y)<1.2&(y-/)<1.5; 


The bolometric IR luminosity, L/r , was derived for all 
IR-detected galaxies, using the most up-to-date models, i.e. 
SED fittin g to U-to-F IR broad-band photometry with Magphys 
Ida Cunha et al.l!T008t) for PACS-detected galaxies, and 24/rm- 
to-bolqmetric_L/R correction using main sequence SED mod¬ 
els by iMagdis et al.l (120121) for MIPS-detected galaxies with no 
PACS data. L/r values were then converted to SFR, to which we 
added a second component accounting for unobscured SFR, in 
order to recover the total SFR (S FR/r vuv)- Assuming S FR/r +uy 
as our benchmark SFR estimate, we derived some relations to 
correct UV and [OII]T3727 luminosities for dust extinction and 
thus obtain a set of consistently calibrated SFR estimators cov¬ 
ering a wide range of redshifts and SFR regimes that will be 
particularly useful in view of the large spectroscopic surveys 
that are currently on-going or will be_carried out in the near 
future (for example, VUDS dLe Fevre et al.ll2015l). VANDELS , 
BigBOSS JSchlegel^etal. 20111) . Euclid (lLaureiis et all 1201 ll) . 
WFIRST (ISnergel et al 1120131) ). 

In Table [4] we summarize all the relations between attenu¬ 
ation and different spectral properties and colours that were de¬ 
rived in the paper. In this respect, the main results of our analysis 
can be summarized as follows: 


1. The UV continuum slope was derived, on spectra for the 
galaxies in the UV sample and from photometry for galaxies 
in the [Oil] sample. The relation between dust attenuation 
(A/rx) and [5 was studied and found to be broadly consis¬ 
tent with literature results at the same redshift, though with a 
larger dispersion with respect to UV-selected samples, which 
we ascribe to our parent sample being mass-selected. Our 
results also suggest that the Calzetti law is valid up to z~3, 
though it tends to over-predict dust attenuation for galaxies 
with the reddest slope. 

2. Using photometric information, the iDaddi et al. (2004!) re¬ 
lation between (B-z) colour and attenuation in galaxies at 


1.4<z<2.5 was updated using A/rx as an independent dust 
attenuation tracer. 

3. In the [Oil] sample (1 <z< 1.6), where the [OII]/13727 emis¬ 
sion line can be detected, we found an anti-correlation be¬ 
tween the rest-frame EW of the line and f3 p hot, which is the 
main result in the paper. Since [i is related to dust attenu¬ 
ation, this relation implies a relation between EWyoi/] and 
dust attenuation for the galaxies in our sample. 

4. We tested the issue of differential attenuation towards, re¬ 
spectively, stellar continuum and nebular emission. Though 
we cannot draw any conclusive statement, our results are in 
line with the traditional prescription of extra attenuation to¬ 
wards nebular lines. 

5. Calibrated SFRs were used to discriminate between differ¬ 
ent SFHs in a quick comparison to SFR estimates from SED 
fitting to UV-to-NIR broad-band photometry (i.e. without in¬ 
cluding FIR information), and the best agreement was found 
with constant star formation, as opposed to exponentially de¬ 
clining. 

6. Finally, we divided our sample into three redshift bins and 
using the calibrations in Table [4] we studied the relation be¬ 
tween SFR and stellar mass. The galaxies in our sample lie 
on a quite close linear relation (cr~03 dex, r Ay ~-0.6) at all 
redshifts, with a slope ~0.7 - 0.8, consistently with other lit¬ 
erature results in the same redshift range. We also find weak 
trends both for the slope and the normalization, with the 
slope becoming steeper and the normalization higher with 
increasing redshift, though the completeness limits and spec¬ 
troscopic requirements of our sample do not allow a robust 
qualification. 
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Appendix A: Derivation of the UV continuum slope 
from spectra 

In Section 15.11 we studied the relation between the UV con¬ 
tinuum slope and dust attenuation. In galaxies dominated by a 
young stellar population, the shape of the UV continuum can 
be fairly accurately approximated by a power law F^A?, where 
F is the observed fl ux ( erg s -1 c rrT 2 A " 1 ) and (3 is the con¬ 
tinuum slope (ICalzetti et al. |1994|) . The wavelength range in 
which B is defined goes from -1250A to -2600A. In this range, 
ICalzetti et alj (119941) defined ten windows to be used to measure 
the continuum, chosen to avoid the strongest spectral features. 
For galaxies at z> 1.6 B was derived through an error-weighted 
OLS (T| 20, by fitting th e av erage fluxes in the ten windows de¬ 
fined by ICalzetti et al.l (1 19941) . in the log(F ,0 - lo g(A) plane. 
More details about the procedure can be found also in lTalia et al.1 
(120121) . We define the slope computed using all the ten win¬ 
dows as Bib- Since not all spectra cover the entire 1250-2600A 
range, we define also two shorter versions of/?: Bis, ranging from 
-1250A to -1800A, and (h<v, - ranging from -1550A to -2600A. 
With respect to the ten original windows, Bis is measured using 
the 7 bluer windows (1st to 7th), while /326s is measured using 
the 6 redder ones (5th to 10th). 

In Table lATl we indicate the number of galaxies for which 
each version of the slope could be computed, depending on the 
wavelength range covered by the spectra. 

In Fig. IA.1I we compare the three definitions of /3 for the 
galaxies whose spectrum covers the entire wavelength range 
1250-2600A. [3 is is always redder (more positive) than Bib, 
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obtain the following relation (r xy ~0.9): 

p 26 = (0.71 ±0.03) x/l 18 + (-0.54 + 0.03) (A.l) 

On the other hand,/Tf, v is systematically slightly bluer (more 
negative) than p 2 6, which is quite intuitive, since /3\& and p 2 (, s 
basically sample the two halves of p 2 (,. In this case, a linear fit 
gives: 


p 26 = (0.86 + 0.04) xp 26s + (0.03 ± 0.06) (A.2) 

with a slightly larger dispersion than in the previous case and 
no dependence on p 2 &. The likely explanation of the discrepancy 
between the three definitions is the presence of a large number 
of closely spaced Fe absorption lines in the 2300 - 2800A range 
dLeitherer et al.ll 19991 : ICalzettill2001 ). For each spectrum in the 
UV sample the appropriate definition of the slope was computed, 
depending on the wavelength coverage. Then all the /?i 8 and p 2 ( )S 
values were scaled to /Te using Eq. A.l and Eq. A.2: these ho¬ 
mogenised values are the ones that were plotted in Fig. [IT] and 
used to derive the relation given in Sec. IP 

As explained in Sec. 15.11 p can be derived also from pho¬ 
tometry. In the [OII]-sample the available photometry forced 
us to adopt a short wavelength baseline and we defined p p hot 
in the range A rest ~ 1500-2600A. As we did with the spectro¬ 
scopic determination of the slope, we used the galaxies in the 
UV-sample to test the effect of different wavelength baselines 
also in the computation of the slope from photometry. We found 
that the values of p p i, ot derived using, respectively, the total base¬ 
line (1250 - 2600A) and a shorter one (1500 - 2600A) are 
broadly consistent with one another, with the short-based val¬ 
ues being on average slightly bluer than the long-based ones, 
though to a lesser extent than in the case of / 3 spec . Since the short 
wavelength baseline is available from photometry for the entire 
galaxy sample, we used directly the values of P p \ w , in the range 
.1 v,'~ ^ 500 2600A to obtain the lelations presented in Sec. IP 

and 15.41 without the need to scale them to the total wavelength 
baseline. 


Fig. A.l. Comparison between p 2 („ Pis, and p 2 (, s , for the galaxies 
whose spectrum covers the wavelength range 1250-2600A. See 
the text for more details about the definitions of ft. Blue dashed 
lines indicate the linear fits to the data. 


and the discrepancy increases for increasing value of the slopes 
(ICalzettill2001l) . 

Table A.l. Number of galaxies for which the different defini¬ 
tions of p could be computed. 


Windows 

1st-10th 

lst-7th 

5th-10th 

Tot. no. galaxies 

35 

33 

46 

only 24pm 

14 

11 

17 

24pm and PACS 

8 

1 

11 


There is a strong correlation between the two definitions of 
the slope. In particular, applying an error-weighted linear fit we 
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